Self Consistent Molecular Field Theory for Packing in Classical Liquids 
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Building on a quasi-chemical formulation of solution theory, this paper proposes a self consistent 
molecular field theory for packing problems in classical liquids, and tests the theoretical predictions 
for the excess chemical potential of the hard sphere fluid. Results are given for the self consistent 
molecular fields obtained, and for the probabilities of occupancy of a molecular observation volume. 
For this system, the excess chemical potential predicted is as accurate as the most accurate prior 
theories, particularly the scaled particle (Percus-Yevick compressibility) theory. It is argued that the 
present approach is particularly simple, and should provide a basis for a molecular-scale description 
of more complex solutions. 
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I. INTRODUCTION 



The disordered packing of molecules at liquid densities 
is aprimary and difficult problem in the theory of liquids 
[HQ]. This problem is typically addressed first by consid- 
eration of model intermolecular interactions of hard-core 
type, interactions that rigidly exclude molecular overlap. 
For those systems, a quantity of primary interest is then 
Boltzmann's available phase space from which follows 
the thermodynamic excess chemical potential discussed 
here. Sophisticated theories, even if esoteric, are avail- 
able mUBQ] f° r tne equation of state of the hard sphere 
fluid. In conjunction with simulation results, adapta- 
tions of those theories provide empirically exact results 
for the hard sphere system [8(. Recent theoretical ac- 
tivity H 0| on the hard sphere fluid emphasizes that 
physical clarity is an important quality of theories that 
might be transplanted to describe more realistic solution 
models. The physical content of available models of pack- 
ing of more realistically shaped molecules is conceptually 
similar to theories of the hard sphere fluid, but the re- 
sultant theories are naturally more complicated than for 
hard spjjeres; Refs. [1 El El El El EE El El El El 
I2f ll2lll22l l2a . l23| give examples of that ongoing activity. 

Recent developments of a quasi-chemical approach to 
the theory of molecular solutions [25| have brought a new 
set of concepts to bear on these problems 26]; these de- 
velopments suggest theories with clear physical content 
and a clear path for progressive improvement. This paper 
pursues these developments further, proposing and test- 
ing a self consistent molecular field theory for molecular 
scale packing. More important than the specific packing 
problem considered here, these self consistent molecu- 
lar field ideas will be carried forward to develop quasi- 
chemical treatments of realistic solutions |25j . 



theory is built upon the relations [25, 2C 
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where p n {R) are probabilities for observing n sphere cen- 
ters in an observation sphere of radius R in the liquid, and 
the K n (R) are well defined equilibrium ratios of concen- 
trations of hard sphere n-complexes with that observa- 
tion sphere with Kq = 1. The quantities K n (R) describe 
occupancy transformations fully involving the solution 
neighborhood of the observation volume. Except in the 
limit of low density, these coefficients are known only 
approximately. Therefore, physically motivated approx- 
imations are required to proceed to specific quantitative 
predictions. 

Our previous study of this problem |26J identified a 
primitive quasi-chemical approximation in which 
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Here v — 4ttR 3 /3 is the volume of the observation sphere, 
fi^ 1 = kT, u(rij) is the interaction between molecules i 
and j (the hard sphere interaction in the present case), 
and £ is a Lagrange multiplier used to achieve consistency 
between the known bulk density, p , and the average den- 
sity in the observation volume. Because of the explicit 
factors of p in Eq. ^ ( will approach the thermodynamic 
excess activity, In £ = (3p, ex with p ex the excess chemical 
potential of Gibbs. The integrals of Eq. [2] are few-body 
inte gral s that can be estimated by Monte Carlo meth- 
ods |26l |. A natural extension of this idea is to approx- 
imate K n (R) on the basis of n- molecule coiffigurational 
integrals that give the low-density limiting quantity, but 
with inclusion of a molecular field PfsCF{f) as 



II. THEORY 



For economy of the presentation, we specifically discuss 
the one component hard sphere fluid. The quasi-chemical 
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inside the observation volume. Averaging of these pro- 
files with respect to the possible occupancies predicts the 
observed density. The consistency sought is then unifor- 
mity of the density, 



6]nK£\R;0tpscF) 
50<PSCf(t) 
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for r inside the observation volume. p is defined by 
Eq. ^ an d in Eq. [S] we have noted that, for hard-core 
solutes, the interaction contribution to the chemical is 

i3^ ex = -hxp Q mm. 




FIG. 2: Example dependence of the density profile on scaled 
molecular field X(3(pscF(r); pd 3 — 0.8. 



FIG. 1: The self consistent molecular field (3tpscF{r) for d- 
diameter hard spheres a spherical observation volume of ra- 
dius d. r/d — is the center of the observation volume, and 
r/d = 1 is the surface. The curves on the bottom panel cor- 
respond, from bottom to top, to reduced densities pd 3 = 0.1, 
. . . , 0.9, in increments of 0.1. The upper panel depicts PfscF 
for pd 3 = 0.9, on a plane through the center of the observation 
sphere. 



We will adopt the convention that the molecular field 
PfSCF(f) be zero at the center of the observation vol- 
ume, and an additive constant be absorbed in the La- 
grange multipliers of the K„ (R; P<psc f ) ■ The molecular 
field PfsCF(f), together with the Lagrange multiplier, 
may be made consistent with the information that the 
prescribed density of the liquid is uniform within the ob- 
servation volume. The density profile for the n-molecule 
case is 1271 



Pn(r) 
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(4) 



FIG. 3: Excess chemical potential of the hard sphere fluid 
as a function of density. The open and filled circles corre- 
spond to the predictions of the primitive quasi-chemical the- 
ory and the present self consistent molecular field theory, re- 
spectively. The solid and dashed lines are the scaled parti- 
cle (Percus-Yevick compressibility) theory and the Carnahan- 
Starling equation of state, respectively. 
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Examples of the results following from these ideas are 
shown in Figs. 11151 These results were obtained from 
a two step iterative procedure from a starting guess 
(3<-Pscf{t) = and the probabilities p n of the primitive 
quasi-chemical theory |26| . With the current approxi- 
mate results, we performed Monte Carlo calculations to 
estimate the densities for each occupancy, and on that 
basis the average density implied by the current field. 
We then updated the molecular field according to 



[Mr)scF] new = [Pv{r)scF] old + / hi 



p(r) 



(G) 



where / is a constant less than one that ensures stable 
convergence of the molecular field; a value of 0.2 was 
found to work here. Convergence is obtained in 20 itera- 
tions of this procedure, typically. Using the field obtained 
above, we then updated the occupancies, reevaluating the 
Kn\R] P^scf) by performing additional few-body sim- 
ulations to evaluate the work associated with turning on 
the molecular field using thermodynamic integration 

K^\R;(3lpscf) = e - J 1 (E"=i^ s cF(r 3 )) A dA 
id o) (i?;0) 

(7) 

where A is a coupling parameter, and (. . indicates av- 
eraging over configurations generated under the influence 
of the molecular field scaled as \(3(Pscf{t*). Using these 

recalculated Kn\R] P^scf), we generated a new set of 
p n , tested for convergence, and so on. This process was 
found to converge within two steps even at the highest 
densities considered. We attribute the observed conver- 
gence to the fact that the starting point, the primitive 
quasi-chemical theory, is accurate for the probable occu- 
pancies. The molecular fields obtained using this method 
were found to converge stably with little difficulty. 

Fig. ^ shows the self consistent molecular fields ob- 
tained using the procedure described above up to fluid 
densities of pd 3 = 0.9, just below the hard sphere freez- 
ing transition. /3ipsCF{r) is a monotonically increasing 
function of radial position from the center of the stencil 
volume to its boundary. This reflects the fact that in 
the absence of the molecular field the hard sphere parti- 
cles tend to build up on the surface of the stencil volume 
to minimize their interactions with the other particles 
(Fig. 0). The molecular field makes the boundary repul- 
sive, depletes the surface density, and homogenizes the 
density within the volume. The magnitude of this repul- 
sive field increases with increasing fluid density. 

The predicted hard sphere chemical potentials as a 
function of density using the primitive and self consis- 
tent molecular field quasi-chemical theories are compared 
to the chemical potential from the Carnahan-Starling 
equation in Fig. |3| The primitive theory works well 
up to pd 3 » 0.35, roughly the critical density for Ar 
and the density region suggested to mark qualitative 
packing changes in the hard sphere fluid [28|: at higher 



densities the primitive quasi-chemical theory systemati- 
cally under-predicts the hard sphere chemical potential. 
The present self consistent molecular field theory sig- 
nificantly improves the agreement with the Carnahan- 
Starling equation over the entire density range. Above 
densities of pd 3 « 0.6, the self consistent molecular field 
theory begins to over-predict the hard sphere chemi- 
cal potential, though the absolute value of the error is 
in marked improvement over the primitive theory. We 
note that the self consistent molecular field theory is 
in closer agreement with the scaled particle (or Percus- 
Yevick compressibility) theory for the chemical poten- 
tial. Fig. 0] shows that the most important deficiencies of 
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FIG. 4: Distributions p n with R—d for densities of pd = 0.35 
(filled circles) and 0.8 (open circles). The dashed lines are the 
primitive quasi-chemical theory of Ref . |2(| , and the solid lines 
correspond to the present SCF theory. Note the marked break- 
away of the n=0 point from the primitive quasi-chemical 
curve, observed before [2^1 . The errors on the high n side 
of these distributions might reflect the fact that the present 
SCF theory doesn't explicitly treat pair correlations. Those 
correlations enter only through the integrals (R; f3ipscF)- 

the primitive quasi-chemical theory are corrected by the 
self-consistent molecular field theory. Note that the self- 
consistent molecular field theory captures the break-away 
at high density of lnpo from the primitive quasi-chemical 
prediction. 

In addition to achieving a uniform density across the 
observation volume, the self consistent molecular field 
also nearly achieves thermodynamic consistency for the 
chemical potential. With the choice of an additive con- 
stant which makes (3^sCF{f) zero in the deepest inte- 
rior of the observation volume, ln£ should approach the 
excess activity of the solvent in the limit of a large ob- 
servation volume. We expect on a physical basis that 
f3(pscF(f) describes an interaction between the interior 
and the exterior of the observation volume across the 
intervening surface. Particularly in the present case of 
short ranged interactions, we expect spatial variation of 
fitpscFir) to be confined to a surface region. Though 
a stencil volume of radius R = d is evidently not large 
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enough to observe that bulk behavior (Fig. , for that 
R = d case we can compare the computed excess chem- 
ical potential with the solvent activity. Fig. El compares 
-lnpo and hiC as determined by the primitive and self 
consistent molecular field quasi-chemical theories. While 
the activity evaluated within the primitive theory signif- 
icantly under-predicts po , with the self consistent molec- 
ular field theory ln£ and -lnpo agree nearly perfectly. At 
the highest densities, there is a slight disparity between 
these two quantities, and the calculated values for ln£ 
are in better agreement with the empirically known f3pf x 
for the hard sphere fluid. 




FIG. 5: Comparison of In £ (with ( the Lagrange multiplier or 
excess activity) against computed excess chemical potential, 
/3p ex — -lnpo, demonstrating the thermodynamic consistency 
of these quasi-chemical theories. The open circles are the 
primitive quasi-chemical theory (Eq. |5J , and the filled circles 
are the present self consistent molecular field theory. 



III. VARIATION WITH CAVITY SIZE 



The related quantity 



4tt P R 2 G{R) = - 



dlnpo 
di? 



(8) 



is of special interest in the theory of the hard sphere fluid, 
and of solubility more generally |^. l2^.l30| . In the present 
quasi-chemical approximation, this is 

4tt P R 2 G(R)^ £ m p m (dlnir ro (#;^scir)ALR) .(9) 

To analyze the required derivative, we consider that the 
radius R is defined in the first place by a bare field [3ipo 
that is zero (0) inside the observation volume and oo out- 
side. Then the full field encountered with the integrals 
Eq.QJis flip — /3ipo + fitpscF- The result now correspond- 
ing to Eq. 0| is 



din K rn (R;f3ipscF) 
dR 



Pm{r;f3<p 



SCF) 



( df3ip{r) 
dR 



d 3 r 



(10) 



Upon separating the contribution from flipo and perform- 
ing the population averaging, this produces the simple 
relation 



4npR 2 G{R) = AttR 2 P - f 

J V 



^|f^V, di) 



With this first population averaging, we emphasize that 
l3ipscF{f) doesn't depend on the occupancy index to. 
The radius derivative d(3ip(r)/dR of the full field can be 
described by a simple formal relation. The relation 



(5p m (r)8p m {r')) 



(12) 



follows from Eq. 21 for each occupancy. Performing the 
population averaging at this stage, we write 



Sp(r 



6/3ip(r') 



- = (5p(r)Sp(r'))= X (r,r') (13) 



and 



8p<p(r) = I x^^OMOdV 



(14) 



Averaging of the functional derivative Eg. 1121 before com- 
posing Eq.^Jis suggestive of the RPA concept of exploit- 
ing an average potential in a linear response function. 
To use Eq. El consider the density change Sp(r') corre- 
sponding to dematerialization of the uniform density in 
a thin shell (R-AR,R). 



dR 



R\ 



X 



(15) 



\f" 



=R- 



where the latter integral is over solid angles covering 
the surface of the ball. We introduce now c(r,r'), the 
Ornstein-Zernike direct correlation function defined by 
X \ r i r ') = S(r — r')/p(r) — c(r,r'). We finally obtain 



G(R) 



1 



/ c(r,r' 

J V 



zR)pd 3 



(16) 



within the present approximation. In the indicated inte- 
gral the r' coordinate is pinned to the sphere surface, and 
r integration is over the interior of the sphere because of 
Eq. The function c(r,r') is the OZ direct correla- 
tion function in the field f3cp including the self consistent 
molecular field, thus for the case of a uniform density 
enclosed in a sphere of radius R with no materal outside. 

It is obvious that Eq. ^3 gives the corrrect answer for 
the case that the solvent atoms have no interaction with 
one another (c(r,r')= 0), and for the same reason this 
formula is obviously correct in the limit of zero density. 
That limiting result gives the second virial coefficient the- 
ory for lnpo- At the initial order in the density 



c(r, r') — exp[— f3u(r, r')] — 1 + O(p). 



(17) 



This relation in the approximate Ea. llbl leads to the cor- 
rect contribution of next order in the density for G(R), 
corresponding the third virial contribution to lnpo- 
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Exact results are also available in the case that the 
observation sphere is sufficiently small, R < f. Then 
po = 1-4tt P R 3 /3, f3ip S cF(r)=lii (l - 4irpR 3 /3) (spatially 
uniform in < r < R, so in the formulation above this 
would be reflected solely in the Lagrange multipliers). 
Direct calculation gives x( r i r ') = M( r — r ') — P 2 : an d 
x - 1 (r,r') = p-^r.r') + 134^73 - Using these re- 
sults in Eg. 1 161 gives the known answer, G(R)= 1 _ i7T 1 pR i 73 ■ 
Tests of other current theories in this regime have been 
given by 0] . 

IV. CONCLUDING DISCUSSION 

The physical content of the present self consistent 
molecular field theory is simple and clear, and this theory 
is as accurate as the most accurate prior theories, partic- 
ularly the scaled particle (Percus-Yevick compressibility) 
theory, for the thermodynamics of the hard sphere fluid. 
The conclusion is that careful attention to the near neigh- 
borhood of a distinguished solute in such a liquid, with 
a self consistent molecular field describing the influence 
of the more distant neighborhood, provides an accurate 
description of packing in dense liquids. Though distinct, 
the hydrostatic linear response theory [To| leads to a simi- 
lar conclusion that good theories of these phenomena can 
be extremely local. The present results address contribu- 
tions essential to quasi-chemical descriptions of solvation 
in more realistic cases , as has been discussed on a con- 
ceptual basis recently |25| . 

The present results provide a definite, and organized 



basis for theoretical study of subsequent solvation phe- 
nomena. For example, consider inclusion of attractive 
interactions between solvent molecule spheres, interac- 
tions secondary to the repulsive interactions. The sim- 
ple estimate c(r) ~ —(3u(r) for distances not too small, 
is consistent with Eq. El But when u(r) at those dis- 
tances describes attractive interactions, Ea. 1 161 then pre- 
dicts that these attractive interactions reduce the mag- 
nitude of G(R). This is a behavior that has been much 
discussed over recent years in the context of theories of 
inert gas solubility in liquid water [2^, |3(J |^| . 

A related but distinct issue is how these packing ques- 
tions are affected by multiphasic behavior of the solu- 
tion, particularly the possibility of drying [3^, |U 0| or 
preferential absorption |3fij in biophysical aqueous appli- 
cations. In such cases, it is attractive to speculate that 
the self consistent molecular field P<fisCF should reflect 
those multiphase possibilities just as it can in pedagogical 
treatments of non-molecular models of phase transitions 
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